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Final Technical Report: Airborne Astronomy Program 
Grant NAG 2-1074 

The Circumxtellar Environment of Low Mass Star Forming Regions 

Harold M. Butner Principal Investigator 

1. Introduction 


Since the 1980's, our understanding about the inter-relationship between the col- 
lapsing cloud envelope and the disk has been greatly altered While the dominant, 
star formation models invoke free fall collapse and a r — 1-- J density profile, other star 
formation models are possible. These models invoke either different cloud starting 
conditions or the mediating effects of magnetic fields to alter the cloud geometry 
during collapse. To test these models, it is necessary to understand the envelope v 
phvsical structure. The discovery of disks, based on millimeter observations around 
young stellar objects, however makes a simple interpretation of the emission com- 
plicated. Depending on the wavelength, the disk or the envelope could dominate 
emission from a star. In addition, the discovery' of planets around other stars has 
made understanding the disks in their own right quite important. Many star forma- 
tion models predict disks should form naturally as the star is forming. In many cases, 
the information we derive about disk properties depends implicitly on the assumed 
envelope properties. How to understand the two components and their interaction 
with each other is a key problem of current, star formation. 


To separate out the two components, it is necessary to observe sources at a variety 
of wavelengths, over a high spatial resolution. As mentioned above, depending on 
the physical conditions in the disk and the envelope, each will contribute dffferrent 
amounts of flux at each measured wavelength. In addition, high spatial resolution 
should allow us to resolve the envelope emission. Envelope emission typically extends 
over many hundreds of AU, whereas the disk will be confined to a much smaller area. 
By combining high resolution multiwavelength data with radiative transfer models, 
one can assess the likely contribution of both the disk and the envelope. The envelope 
contribution dominates the far-infrared emission, with as much as 80% of the total liux 
coming from the cool dust found in the outer portions of the envelope. This emission 
can be detected and mapped with high resolution observations, such as the Yerlo.-s 
Far-InfraRed (FIR) Camera on board the Kuiper Airborne Observatory! KA0 1 In 
1993 and 1995, we obtained far-infrared (60. 100, 160, 200 micron) images for a 
number of young stellar objects (YSOs) with the KAO prior to the KAO s shutdown. 


In the proposal covered by the grant (NAC-2-1074). we proposed to analyze the 
data, acquire additional data at other wavelengths for models, and to develope source 
models for the YSOs under study. Based on the results of these models, we hope to ^ 

define the range of physical conditions that are expected in young stellar systems ofj^ (j 
different ages and masses. Through the Airborne Astronomy program, we were able ' 
to accomplish man} - of our project goals. MAY 1 0 1 
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II. Observational Data Acquired for Source Models 

Far-Infrared Images: In FY 93 and FY95, we were awarded two flights each 

year 1o study the envelope arounds low mass protostars in Taurus. Working with 
the Ycrkcs KIR. camera, wc obtained data on a total of 13 sources. The Ycrkes FIR. 
camera, contained an t>(.) element, array. Our standard observing method was to dither 
the array around the desired source position, and to take a series of images for later 
co a.d ding. Initial data reduction indicated that most of the sources were compact., as 
would be expected for low luminosity protostars surrounded by an irifalling envelope. 
Such objects will not tend to heat up a lot. of the envelope, and what dust is heated 
will be concentrated towards the center. However, two sources (1.1527 and L155J.NE) 
appeared to show extended, emission. This was particularity interesting since these 
objects may be the the youngest objects in the sample. The claim of youth is based 
on their spectra] energy distributions (SF.D.s), which are. peaked very strongly in the. 
FIR. Such a peak implies that the object is still deeply enshrouded in dust. 

Our initial data reduction was based on smoothing the data with a relatively 
broad smoothing function (i.c. a Gaussian beam). The smoothing resulted in the 
data being far from the diffraction limit. To address this situation, we undertook 
[ with the Yerkcs group) to re-reduce the data using improved convolution methods. 
Due to the desire, we had to spend considerable time to achieve the diffraction limit. 
'I 'hose new results are still being analyzed, however, they confirm the detection of 
extended emission around L1551NE and L1527. The final models of the KAO results 
are still in preparation. 

To supplement the KAO data, we obtained observations at related wavelengths 
for a number of the sources in our KAO dataset. 

Mid-Infrared spectra: A large sample of T Tauri stars was studied by Cohen 

and Wittebom in 1965 at 10 pm. To extend this survey, we obtained R=200 spectra 
for several T l'auri stars, and a more extended sample of Herbig Ae/Be stars. These 
stars hate a variety of 9.8 //m silicate emission. This emission feature reveals that 
optically thin silicate dust is in the system. Such dust is similar (spectroscopically) 
to dust seen around some Vega-like objects like Beta Pie. 

The passage of cornel Hale Bopp provided an opportunity to measure the spec- 
troscopic properties of comets for comparison with our T Tauri arid Herbig Ae/Be 
sample. The comet spectra was found to have some features that were time variable 
The comet, spectra reveal the presence of grain componenents that might l>e pre-solar 
(based on comparison with IDP spectra). We hope to tie the Comet Hale Bopp work 
to our FIR/ MIR studies m our final mode/s. 

Mid-Infrared images: Using the MIRAC2 (Mid-Infrared Array Camera), we 

imaged several of the embedded YSOs at 10 i. These observations provide high 
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sensitivity high spatial observations of the circumstellar environment on size scale* 
of 1 SO AU. Most of the observations arc consistent with the mid-infrared emission as 
coming from a disk. Extended emission is possible, since we found several examples 
in related studies of Herbig Ae/Bc stars. However, in general, the T Tauri / embedded 
low mass stars are too cool to heat an extended envelope unless some unusual grain 
heating mechanism is invoked. 

We obtained 10 and 20 /im images of the Trapezium cluster, with its proplyds, 
which we arc currently analyzing to provide a. comparison with our isolated stars. 
Initial results are being published in Deutscb ct. al (1999). 

Submillimeter/millimeter photometry: We focused on Taurus embedded 

YSOS and the Vega-like stars. In the case of the Taurus sources, which includes our 
KAO sample, we have obtained a complete survey of the. 450-2000 tun emission using 
the JCMT 15 motor telescope. For the Vega-like stars, wc used the CSC 10m and 
IRANI 30m telescopes to get 300 and 1300 /;m photometry. The ‘data will be written 
up in two papers. 

Submillimeter/Millimeter Interferometry: 

A number of our sample stars have data already available in the literature. For 
t hose sources not already observed, we obtained 1.3 mm interferometry data for the 
brightest 16 YSOs in our Taurus sample. These observations extended a. previous 
survey we had done. We detected 14/16 sources (Wannior ei al. 1999 in prep). Of 
particular interest, we discovered evidence for a couple of binaries in our sample 
(‘L.1551NEK This dataset is to be combined with published interferometer 2.7 .mm 
data reported in several surveys. Wc also were awarded time to study a sample of 
Taurus YSO$ with the JCMT-CSO Single-Bearn-Interferomotcr, and got data on DC 
Ta u , a m om$ ot h or sources . 

Summitry: Wc have obtained the complete $ED from 10 microns out to 

1.3mm for all of our sources. We have the FIR imaging data, processed to reveal the 
maximum angular resolution possible, which allows us to model the disk lb model 
the disk, we have the high resolution millimeter interferometry data. 


L AV 

inA 


Klf 1 

11 u. 


tnnnotmn a 
JCUOCOOOOH 


1 A 

14 


P r\ 4 

. U4 


1 


111. Modeling Techniques Developed for Source Models 


At the start of our KAO research, we had one primary modeling technique; a 1-D 
radiative transfer code that allowed us to explore a variety of density gradients and 
grain properties. We also developed a way to including disk emission in our model. 
With these models, we ran explore many types of disk/cnvelope models. 

We identified two other radiative transfer codes that would be useful for our 
models: 

]..) Inclusion of Polycyclic Aromatic Hydrocarbons (PAHs) and small grain heat- 
ing for the models. This code, developed by Siebenmorgcn, Krugel, and Mathis 
(1992), allows us to do spatial models of the PAII emission and small grain emission. 
Such effects are critical to understanding the emission from a number of the Herbig 
Ae/Be stars we have under study. 

2.) A fast 2-D code: Developed at Jena by Dr. Men shchikov, it can handle 
both deeply embedded sources as well as visible disks. It has been applied by Dr. 
MerPshchikov to LI 55.1 and 11 JL Tan. We intend to apply it to our more deeply 
embedded sources (Li 551 N E and L 1 52 < ) . 

Overall, there are three important parameters we examine: 

1) The envelope geometry and temperature structure. 

This is best done with the SED and spatially resolved data at far-irifrared/snlv 
; n i 1 1 i meter wa vel engths. 

2) The disk geometry and temperature structure. 

This is most sensitive to the mid-infrared emission, due to the role of small grains 
or PAHs. High resolution data at 10 microns provides the best constraints. 

3) The clust properties in both the envelope and the disk. 

The models allow a tradeoff (to a limited extent) between dust properties and 
the inner stucture of the disk /envelope system. However, the presence or absence of 
silicates is strongly limited by the 10 micron spectra. 
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IV'. Results to (late 


Combined with our iuitia.1 K AO program, we find: 

1. ) The vast, majority of embedded YSOs in Taurus are compact at 100 microns. 
The models most consistent with our data and other observations are either dominated 
by disk emission (for sources like GO Tau) or envelopes that have relatively sleep 
density gradients ( r _1 -5to-2.0^ 

2. ) The submillimet.er/millimeter photometer suggests that envelope models 
(such as those of Hartmann ft at. (1903) are very successful. However, in a few 
cases) they are clearly missing a component, in their models. Disk emission plays an 
important role and must, be considered when predicting the overall emission. I Butner 
ct ul. 199 9b (in preparation)). 

o) Tn t.hc two cases where we seem to have extended emission, wo have to investi- 
gate other possible source models than a Shu collapse. This work is still in progress. 
We will considering models such as those of Hartmann cl ul. t 1995 l , which invoke a 
flat sheet, rather than spherical collapse. 

V. Future Work 


We still cloim\ the final modeling for the revised FIR maps, and anticipate having 
two or three papers done by the end of the year. These papers will present our data 
and ihe complete se" of parameter spa.ee that we explore. Tn particular, we. wish 
to assess the role of flat sheet collapse models (see point. 3 above) in explaining our 
Taurus FIR data. 
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VI. Papers in Preparation or Published 


The following papers were supported in part by the Airborne Astronomy Grant 
NAG 2-1074, The support provided salary, travel funds, or equipment used for data 
reduction and analysis. 

Refereed Articles 

H. J. Walker, V. Tsikoudi; C. A. Clayton, T. Gcballc, D. H. Wooden, k H. M. Butner, 
“The Nature of the unusual source IRAS 18530+0817/ 1997. A&A, 323, 412. 

H, M. Butner, H. J. Walker, D. II. Wooden, k F. C. Wilteborn, *1997/’ Examples of 
Comet- Like Spectra Among j3 Pic-Like Stars, in Astronomical and Biochemical 
Origins and the Search for Lift in the Universe, ed. C. B. Cosmovici, S. Bovvyer, 
k D. Wcrthimer ( Bologna: Editrice Compositori), 

C. H. Moriarly-Schieven k H. M. Butner, “A Sub-Millimeter- Wave 4 Flare ? From GG 

Tan?, 7 ' 1997, Ap, 7., 474, 768. 

Selected Unrefereed Articles 

A. Natta, k II. Butner, * Resolving Disks in YSOsV 1996, in Infrared Space Inter - 
fcromttvy , ed. A. Alberdi, T. do Graauw, C. Eiroa, C. J. Schalinski, and H. 
Tli ro nso n f S p ace S deuce Re v i ews ) , in p ross . 

Papers in Preparation 

H M. Butner. (J. II. Mori&rty-Schievcn. P. Wannier, “Spectral Energy Distributions 
in Taurus/* 1999, Ap.J.. in preparation. 

H. M. Butner, G. M. Moriarly-Schieven, ei a l. . “FIR Imaging of Low Mass Sources 
in Taurus/ 1999, Ap.J,, in preparation. 

L. K. Deutsch, M. Kassis, N. Smith, J. I. Hora, G. G. Fa^io. H. M. Butner, W. F. 
Hoffman. A. Dayal, “Wide-Ficld Thermal Imaging of the Orion Nebula at High 
Spatial Resolution/’ ] 999, Ap.J.< in preparation. 

P. Q. Wannier, J. Keene. G. Moriarty-Schieven, and H. M. Butner. *\l..3 MM inter- 
ferometry of Taurus Sources/’ 1999, Ap. in preparation. 

D. II. Wooden, H. M. Butner, k K. (\ Witteborn, “A Mid-Infrared Survey of T Tauri 

Stars.” 1999, Ap. 7., in preparation. 
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